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I SUMMARY 

i 

cTDnM^MrtpTM me Presented which allows the modeling of a moving boundary with 
STRAN NOLIN cards. Various aspects and limitations of the approach are 
plained. Recommendations are given as to the procedure to be used in imDle- 
nting the method. K 

i 


| INTRODUCTION 

I 

Situations occasionally arise when the boundary conditions of a theoretical 
pel change as a result of the response of the model. The case under consider- 
[ion is that where the boundaries of the model move as a function of rigid body 
s placement. The problem which motivated the development of the technique which 
, are presenting is that of a shaft under both axial and transverse loading movinq 
jtween two stationary bearings. Figure 1 is a diagram of such a shaft. 

The reactions at the bearings are dependent on the applied loadings and the 
brtial forces generated by transverse motions of the shaft. The applied loadings 
F * nowri beforehand. The inertial loadings of the shaft are unknown since the 
ansverse motion of the shaft is an unknown. Therefore, the reactions of the 
arings cannot be predetermined before the start of the analysis. 

Constraints cannot be directly applied to the model to account for the 
arings since such constraints would be useful for only a short time of the 
alysis covering a small axial displacement of the shaft. It is possible to solve 
s problem by dividing the total time of the analysis up into many small time 
:rements each with a different constraint condition. The idea here would be to 
art the analysis with the bearings in one position; run the problem for a short 
Tie; stop the analysis; change the constraints; and restart the problem with 
itial conditions equal to the state of the model at the end of the first time 
Spent. This process is repeated until the shaft has moved as far as desired 
is technique, while being possible, is awkward from the standpoint of data 
ndling, the number of runs required, and the total number of da vs needed to 
tain an answer. 


SYMBOLS 

distance between grid points 

* F 2 nonlinear loads at grid point accounting for moving boundary 
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applied time dependent axial and transverse loads 
index number of grid point 

l 

maximum slope of nonlinear loading curve 

i 

stiffness of bearing 
scale factor 

grid point axial and transverse displacements 

axial and transverse coordinates of model 
integration variable in axial direction 

DISCUSSION j 

The approach which accounts for the moving bearings treated the bearingsj 
as stiff springs which moved along the model of the shaft as the shaft moved wiij 
respect to the bearings. These springs were fabricated artificially by the use 
of the nonlinear load option in NASTRAN. Each grid point along the shaft model 
over which the bearings could pass had two nonlinear loadings applied to it for 
each bearing. Figure 2 shows a diagram of these loadings for a typical grid poi 
on the shaft model. For simplicity we consider the case of only one bearing. 

The force F 2 is shown as a positive ramp function with a slope equal to 1 

bearing stiffness and is a function of the displacement u x of the shaft. The 

second force F-j is shown as a negative ramp function with a slope equal to the 

bearing stiffness and is a function of the shaft displacement u x added to the 

transverse deflection of the grid point u y . Note here that since u y is several 

orders of magnitude smaller than u x , the magnitude of u x + u y is mainly determii 

by u alone. Both F-, and F £ are zero until u x is greater than x^-,, the distan 

from the bearing at the start of motion to grid point i-1. F-j and F 2 are equal 

and opposite and, hence, negate each other when u x is greater than x i+1 , the 

distance from the bearing at the start of the motion grid point i+1 . Therefore 
the loadings F 1 and F 2 add to zero at all times except those times when the 

bearing is between grid points i-1 and i+1. 

Now consider a case when the bearing is between grid points i-1 and i+1. 
Figure 3 illustrates such a condition. 

NET LOAD = F-| + f 2 = _K ^ U x + u y ' x i-l^ + K ^ U x " x i-l^ 

NET LOAD = -Ku y 
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s can be seen, the net load on the grid point i is equal to -Ku which is a 

•estoring force proportional to the bearing spring constant and the lateral 
leflection of the shaft at i. Thus, the combination of two nonlinear loadings at 
results in simulating a bearing spring which turns itself on when the bearing 
s in the vicinity of grid point i and turns itself off when the bearing is out- 
side the vicinity of grid point i. 


It is possible to allow the spring constant of the bearing to vary as a 
unction of its distance from a grid point. There is some advantage to this 
;ince when a bearing is between two grid points, its stiffness ought to be shared 
n some way by the two grid points. See figure 4. 


We have used the proportional sharing scheme where the stiffness of the 
gearing is distributed in proportion to that of a load shared between two ends of 
i pinned ended beam. To establish the loading curve for this kind of stiffness 
sharing scheme, notice that 
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(2a - x) 2 

(4) F - Kj|a 2a J 
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Equations (3) and (4) establish the shape of the loading functions for ! 
the proportional stiffness loading scheme. Other stiffness sharing schemes 
could be devised and their loading curves easily arrived at in a similar manner! 
to that shown here. Figure 4 shows a plot of the nonlinear load function for 
the proportional stiffness sharing scheme we used. 


The scheme depends completely on accurately determining the difference 
between the quantities K(u x + u ) and Ku x - Theoretically, this presents no 

difficulty. If the value of u x is six or seven orders of magnitude larger than 
u , a numerical difficulty arises since a computer may only work to eight signi- 
ficant figures using single precision arithmetic. Even if double precision 
arithmetic is used, there will come a point where large values of u x start to 

degrade the numerical accuracy of the differencing operation K(u x + u y ) - Ku x » 


There is a way around this numerical difficulty. The objective of the 
differencing operation is to obtain the quantity Ku y . If the u x quantity is 

multiplied by a scalar constant the differencing operation remains unchanged. 
K(su x + u y ) - Ksu x = Ku y 

Here u was multiplied by the scalar s. To remove the numerical difficulty, 
the scalar constant is assigned such that u x is no more than three or four 

orders of magnitude greater than values of u y which can possibly be of interest 

su cannot be too small or else the larger values of u y will start to affect 

the first few significant figures of the quantity su x + u y . This would result 

in inaccuracies in the determination of the spring constant since the spring 
constant should be a function of u x only and not su x + u y . Thus, some judgment 

must be used in the determination of s; that is to say it cannot be made some 
very small number without considering the values which u y are expected to be lr 
the neighborhood of a bearing. 
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When using the scalar multiplier the nonlinear loading diagram changes 
jch that the nonlinear loads start to act at sx^ rather than at s.. ^ since 

ie lock up function or independent variable is su x rather than u x - Figure 5 

lows the forcing function set up using the scalar multiplier. 


The equation for the nonlinear load as a function of su or su + u 


ecomes 
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The slope of the nonlinear force curve with respect to the scaled 
coordinate remains the same as the slope of the force curve using an unsealed 
fcoordinate. 

f Figure 6 shows a typical result for the transverse displacement of a grid 

point in the analysis of a shaft moving between two bearings. The results show 
the expected suppression of transverse displacement while the bearings are in the 
vicinity of the grid point. 

The application of nonlinear loads raises questions as to the numerical 
stability of the solution. The numerical stability will depend on the time step 
used and the maximum slope of the nonlinear force curves. If the moving boundary 
is a rigid one, some compromise must be made on the stiffness of the moving 
boundary. An infinite stiffness will always result in an unstable condition. 

In our case, we ran several short problems with our model using what we considered 
as a reasonable time step. We reduced the slope of the nonlinear load curves 
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until we obtained a stable solution and then reduced the time step until the 
solution appeared reasonably converged. This trial and error approach may not 
always work out, but it did in our case. We were able to arrive at an acceptabl< 
time step and loading slope within a few days. It is best when setting up the 
tables for the loading curves to set the maximum slope in the table to unity and 
to use the NOUN cards to scale the slope. Thus, the maximum slope values will 
appear directly on the NOUN cards. The tables are tedious to construct as then 
is a table for each grid point on the model over which the moving boundary must 
pass, whereas it is relatively easy to change the slope values on the NOUN cards 

If the time steps are large, the slope of the nonlinear load curve must be 
very small in order to give good convergence and numerical stability. In such a 
case, the moving boundary cannot be accurately modeled using this method. This 
is unfortunate but as far as we can see, there is no alternative for solving i 
the problems of convergence and stability other than to reduce the time step 
or the loading slope. 
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Net Load = F] + F 2 * -Ks(u x +u y -x f l ) +Ns(u J( -x^_ 1 ) = -Ksu 

FIGURE 3. - NONLINEAR LOADS WITH BEARING NEAR 
GRID POINT 











